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ABSTRACT 


A fluid  mechanical  model  is  developed  to  assess  the  performance 

of  a rocket  that  is  propelled  by  the  absorption  of  radiant  energy  from  a 

remotely  stationed,  repetitively  pulsed  laser.  The  model  describes  the 

flow  within  a conical  nozzle  that  is  subjected  to  point  energy  depositions 

at  the  apex  of  the  cone.  A similarity  solution  is  obtained  and  the  specific 

impulse  and  energy  efficiencies  that  may  be  achieved  with  such  a device 

are  determined.  Fluid  mechanical  constraints  limit  the  range  of  pulse 

repetition  rates  that  may  be  utilized.  Preliminary  design  considerations 

indicate  that  a specific  impulse  of  800  seconds  or  greater  may  be  achieved 

with  both  a laboratory  and  a full  scale  device.  A two  pound  laboratory 

rocket  can  be  accelerated  at  10  g's  with  a 15  joule  laser  pulsed  25,000 

times  per  second.  A one  ton  rocket  will  require  a megajoule  laser 

operating  at  350  pulses  per  second  to  achieve  an  eq  livalent  acceleration. 

A laboratory  experiment  to  test  the  theoretical  model  using  multiple 

n 

CO2  TEA  lasers  is  also  designed,  and  a test  plan  to  compare  theory 


with  experiment  is  outlined. 
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I.  INTRODUCTION 


In  recent  years,  several  authors  have  discusser!  and  analyzed 
the  possibility  of  beamed  laser  energy  for  rocket  propulsion,  often  with 
specific  reference  to  the  application  of  high  power,  ground  - based  lasers. 
The  concept  is  deceptively  simple:  provide  a high  energy  density  for 
propulsion  without  the  encomberance  of  a massive  on-board  power  supply 
by  absorbing  radiation  from  a remotely  stationed  high-power  laser. 

Since  the  radiation  absorbing  propellant  may  be  a high  temperature  plasma, 
the  specific  impulse  can  be  very  large.  The  achievable  thrust  may  only 
be  limited  by  the  available  laser  power,  and  with  a remote  energy  source 
large  payload /vehicle  weight  ratios  are  possible. 

A series  of  experiments  to  determine  the  specific  impulse  and 
thrust/laser  power  that  can  be  obtained  with  existing  laser  systems  are 
described  in  Ref.  4.  Steady-state  simulation  experiments  were  performed 
in  a vacuum  chamber  with  solid  propellants,  and  pulsed  laser  propulsion 
along  with  the  laser  - powered  pulse  jet  concept  is  introduced.  A steady- 
state  or  CW  laser  propulsion  system  is  a system  whose  thrusf  remains 
constant  in  time  while  the  laser  beam  continuously  provides  the 
energy  source  for  converting  propellant  mass  to  exhaust  kinetic 
energy.  It  was  found  in  Ref.  4 that  a high  ratio  of  thrust  to  laser  power 
can  be  obtained  by  simply  using  the  laser  to  vaporize  a solid  surface. 
However,  in  order  to  obtain  high  specific  impulse  it  is  necessary  to  add 
energy  to  the  vapor  in  a stable  manner.  The  heating  of  a gas  by  external 
radiation  downstream  of  a nozzle  throat  was  found  to  be  inherently  unstable 
when  the  gas  is  initially  weakly  ionized  and  absorbs  radiation  via  inverse 
Bremsstrahlung.  The  stability  of  laser -heated  flows  both  upstream  and 


downstream  of  nozzle  throat  is  not  adequately  understood  and  is  a 

7 

very  complex  issue.  However,  it  appears  that  stable  heating  of  a 
propellant  in  a steady  state  manner  may  best  be  accomplished  by  heating 
the  gas  upstream  of  a throat  such  that  the  beam  direction  and  the 
propellant  flow  direction  are  the  same.  This  would  require  a laser 
window  in  the  absorption  chamber  that  will  tolerate  transmission  of 
significant  laser  intensities  along  with  high  pressures  for  long  periods  of 
time. 

The  alternative  approach  to  laser  propulsion  that  circumvents  the 

stability  problem  is  to  utiliz.e  a pulsed  laser.  The  techniques  for  obtaining 

large  thrust  and  specific  impulse  with  a pulsed  laser  are  an  outgrowth  of 

8-12 

various  experimental  and  theoretical  programs  in  laser  effects.  When 

a high  power  pulsed  laser  is  focused  to  a high  intensity  in  a gas  or  on  a 
solid  surface,  a high  temperature,  high  pressure  plasma,  which 
propagates  up  the  laser  beam,  is  initiated.  Provided  the  pulse  is  suffi- 
cieitnly  short  that  the  high  pressure  gas  remains  in  the  vicinity  of  a surface 
or  nozz.le  wall,  this  method  may  be  an  efficient  propulsion  mechanism. 

The  propulsion  system  operates  in  a way  similar  to  detonation  propulsion 

systems  that  have  been  proposed  for  use  in  high  pressure  environ- 
1 3 

ments.  Periodic  'explosions'1  in  the  nozzle  transfer  the  detonation 
or  laser  energy  to  the  working  fluid.  The  two  most  significant  potential 
advantages  afforded  by  a pulsed  laser  propulsion  system  over  a CW 
laser  propulsion  system  are  1)  simplicity  in  engine  design  as  a result  of 
permitting  the  laser  beam  to  enter  the  nozzle  via  the  exhaust  plane,  and 
2)  elimination  of  constraints  resulting  from  plasma  instability.  However, 
the  power  conversion  efficiency  (efficiency  of  converting  laser  power  to 
power  in  the  rocket  exhaust)  must  be  determined.  In  Ref.  4 a low 
power  conversion  efficiency  was  obtained  because  the  pulse  time  of 
the  laser  was  too  long.  In  addition,  with  pulsed  laser  propulsion. 
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thrust  is  obtained  when  laser  energy  is  converted  to  kinetic  energy  by  a 
continuously  weakening  shock  wave  traversing  the  propellant  gas.  The 
relative  efficiency  of  generating  thrust  in  this  manner  is  not  known  a priori 
to  be  the  same  as  when  converting  laser  power  to  thrust  in  a steady  process. 

The  purpose  of  the  present  report  is  a detailed  study  of  the  fluid 

mechanics  of  pulsed  laser  propulsion.  As  a result  of  this  theoretical  study 

the  laser  requirements  for  an  experimental  test  of  pulsed  laser  propulsion 

concepts  are  specified,  and  a suggested  experiment  is  presented.  The 

objective  is  to  determine  the  relative  efficiency  of  a pulsed  laser  propulsion 

system  compared  to  a CW  system  and  to  calculate  the  specific  impulse  and 

thrust  as  a function  of  laser  power,  pulse  repetition  frequency,  ambient 

conditions  and  propellant  mass  flow.  The  nozzle  configuration  is  taken  to 

be  an  idealized  extension  of  the  concept  introduced  in  Refs.  2 and  4.  A 

2,  4 

schematic  of  the  single  pulse  nozzle  configuration  is  presented  in  Fig.  la. 
The  nozzle  walls  focus  the  incoming  beam  to  yield  a breakdown  in  the  air  at 
the  focus.  With  a short  laser  pulse,  the  resulting  shock  becomes  a blast 
wave  which  propagates  to  the  nozzle  exit  plane,  converting  all  of  the  high 
pressure  gas  behind  it  into  a force  on  the  nozzle  wall.  This  nozzle  was 
designed  for  single  pulse  operation  only.  Therefore,  no  considerations  of 
propellant  supply  were  necessary.  In  this  report  the  fluid  mechanics  of  a 
repetitively  pulse  laser  propulsion  system  is  analyzed,  and  thus,  the  fluid 
dynamics  of  the  propellant  feed  system  is  included.  The  configuration  to 
be  analyzed  is  shown  in  Fig.  lb.  The  nozzle  drawn  with  a solid  line  is  the 
parabolic  self  focusing  nozzle.  However,  for  simplicity  this  nozzle  is 
replaced  by  a conical  nozzle  which  is  shown  dashed  in  the  figure.  The 
angle  of  the  cone  is  chosen  such  that  the  exhaust  gases  leave  the  exit  plane 
at  the  same  angle  relative  to  the  thrust  axis  as  with  the  parabolic  nozzle. 

The  beam  is  assumed  to  be  focused  externally  so  that  the  focusing  angle 
equals  the  cone  angle.  The  propellant  is  treated  as  a steady  source  flow 
entering  at  the  apex  or  "throat"  of  the  conical  nozzle,  and  periodically 


Aluminum  Parabolic  Shell  of  Revolution 

\ 


(b) 

Fig.  1 a)  Single  Pulse  Propulsion  Concept  as  Introduced  in 
Ref.  4 

b)  Multiple  Pulse  Propulsion  Concept  and  Equivalent 
Conical  Nozzle 


laser  induced  blast  waves  are  ignited  at  r = 0 where  r is  measured 
from  the  apex.  In  Sec.  II  the  theory  for  single  pulse  propulsion,  where 
the  propellant  gas  expands  into  a vacuum,  is  presented.  Sections  III 
and  IV  describe  the  extension  of  the  analysis  to  multiple  pulses  to  determine 
the  specific  impulse  and  thrust  as  a function  of  pulse  repetition  frequency. 

A discussion  of  operational  limitations  such  as  the  effects  of  the  cone 
angle,  a finite  exit  plane  back  pressure,  propellant  plenum  chamber 
conditions  and  nozzle  curvature  is  included  in  subsequent  sections.  The 
design  of  an  experiment  to  test  the  theory  of  pulsed  laser  propulsion  is 
presented  in  Section  VII. 
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II.  SINGLE  PULSE  THEORY 


Consider  the  steady  propellant  flow  through  a conical  nozzle  of 
solid  angle  O.  Utilizing  spherical  symmetry  within  that  solid  angle,  we 
obtain  the  gas  density 


P 


1 


* „ 2 

p (r  ) / 


2 

r 


(i) 


where  p is  the  density  at  sonic  conditions,  r is  the  spherical  radius  and 
r is  the  radius  of  the  hypothetical  "source"  of  strength  M, 


* 2 

M = p u (r  ) Q , 

where  u is  the  sonic  gas  velocity.  The  gas  velocity  in  the  supersonic 
portion  of  the  nozzle  is  assumed  to  be  the  limiting  velocity  of  the  gas  at 
zero  temperature. 


Equations  (1)  and  (2)  represent  the  gas  density  and  velocity  downstream  of 

the  source  radius  r . There  is  a local  failure  of  the  equations  for  r sO(r  ) 

However,  mass  flux  is  conserved  within  this  radius. 

At  t = 0 we  shall  (with  a pulsed  laser)  deposit  energy  E at  r = 0. 

This  energy  will  generate  a shock  wave  which  propagates  with  velocity  Vs  (t) 

spherically  outward  to  radius  Rg  (t).  Equations  (l)and  (2)  represent  the 

conditions  upstream  of  the  shock.  Therefore,  the  shock  propagates  into  a 

2 

gas  whose  density  ~ 1/r  . It  has  already  been  assumed  that  u is  much 

X 

greater  than  the  upstream  speed  of  sound.  Hence,  if  we  assume  that 
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V >:>  u the  shock  may  be  considered  strong  in  the  traditional  sense  of 

14 

the  strong  blast  wave.  The  density,  pressure,  and  gas  velocity  imme- 
diately behind  the  shock  become 


p'  (r")  /R 


(t)  , 


PS(U 


_J_  hZL  * 

(Y  + Y + 1 P 


( r ')  V 2 (t)  /R  2 (t)  , 
s s 


and 


us  (t) 


vs  (t)  , 


respectively. 

In  order  to  obtain  the  density,  pressure  and  velocity  everywhere 
between  r - 0 and  r R (t),  we  must  solve  the  inviscid  Euler  equations. 
The  conservation  of  mass,  momentum  and  energy  are,  respectively: 

& P , & (p  »)  , '!■  p u 

&t  Jr  r ’ 


j u 

a t 


a u i 

a r p b r 


0 , 


(3) 


and 


sip/.’)  . u a (p/o') 

at  b r 


We  seek  a solution  to  Eqs.  (3)  for  the  blast  wave  propagating  into 
nonuniform  gas  that  is  similar  in  r/R  (t).  That  is,  we  assume 


p ps(t)  f(’l), 

p ps (t)  gCn), 


a 


-8. 


and 


where 


u = u (t)  h(T|)  , 
s 


■H  = r/R  (t)  . 
s 


Substituting  the  self- similar  profiles  into  the  Euler  equations,  we 
determine  the  shock  motions  for  which  the  variables  "H  and  t may  be 
separated: 

R (t)  = Atm  , (4) 

s 

where  A and  m are  arbitrary  constants.  The  thermal  and  kinetic  energy 

in  the  shocked  gas  must  be  equal  to  the  energy  deposited  at  r = 0. 

R (t) 
s 

E =/  (CyT  + iu2  ) P n r2  dr  (5) 

o 

From  Eqs.  (4)  and  (5),  it  follows  that 


m = 2/3 


and 


A = E/ 


[rV 


fi  I, 


1/3 


whe  re 


1 

8 / (g  + f hZ  ) q2  d r] 


9 ( y + 1)  V(y+1)  V (Y  - 1) 


Equation  (4)  is  valid  for  the  time  period  over  which  the  shock  is 

strong.  The  limiting  time  is  the  time  at  which  the  shock  velocity  decays 

to  the  upstream  gas  velocity.  This  is  denoted  by  t and  is  determined 

max 

from 


V 

s 


(t  ) = u 
max  i 
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or 


3/2 


r 


t 

max 


8 E 


(x_L\ 

\y  + l ) 


27  n I,  P (r  ) (u  ) 


(6) 


Having  successfully  separated  the  T)  and  t variables,  the  resulting 
equations  for  the  conservation  of  mass,  momentum  and  energy  become, 
respectively. 


t / V f 


n f 


+ (y  + i)  - 


2 h hf 


2 \ Y - 1 


1 '^L\hf  L _ 


Y - 1 


(Y+  1) 


(7) 


and 


(y  - 3/2)  (Y  + 1)  g + (h  - ^ ^ TL  i ig  . f 


&±\  0 . 


Equations  (7)  are  integrated  from  T)  = 1 to  r\  = 0,  subject  to  the  condition 
that 

h (1)  = g (1)  = f (1)  = 1 . 


The  density,  pressure  and  velocity  profiles  are  illustrated  in  Figs.  2,  3 
and  4 respectively.  These  solutions  will  be  utilized  to  determine  the 
performance  of  a pulsed  laser  propelled  rocket. 


Density  Profile  Behind  Strong  Shock 


I 


vtt.1  * .'■iPl^ii3JiiL||iiii..i» 


2- 


7=1.67 


-r/Rc(t ) 


I 


Before  extending  the  fluid  mechanical  model  to  include  pulse 


sequencing,  we  wall  first  define  and  relate  the  time  average  quantities  in 
a notation  similar  to  that  used  for  conventional  rockets. 

Consider  a pulsed  laser  propulsion  system  in  which  mass  M is 
released  per  cycle  and  t is  the  time  between  pulses  (in  contrast  to  the 
pulse  duration).  The  time  average  mass  flux,  thrust  and  power  are 
defined  by 

M M/t  . 


f / 


P 


T tt 


and 


P 

f / 


P dt  , 


respectively.  The  specific  impulse,  I , is  defined  in  terms  of  the  time 
average  thrust  and  mass  flux 


I 


sp 


M g 


where  g is  the  acceleration  of  gravity.  In  addition  to  the  time  average 
quantities,  it  is  necessary  to  define  a mass  average  quantity.  Denoting  the 
nozzle  exit  plane  velocity  of  the  mass  element  dm  as  t^,  we  express  the 
mass  average  exit  plane  velocity  as 


M 


Ue  M f 


o 

-15- 


u dm  . 
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(8) 
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We  seek  to  determine  the  optimum  thrust  T for  the  minimum  time 

average  laser  power  P.  In  the  limit  of  zero  temperature  at  the  exit  plane, 

the  instantaneous  thrust  T is  related  to  u by 

e 


and  the  instantaneous  laser 


o 

power  Pis 


M u , 
e 

related  to  the  kinetic  energy  by 


From  the  above  relations,  the  ratio  of  thrust  to  power  is  expressed  as 


t 


and  the  mass  average  exit  plane 
via 

I 

sp  e 

Eliminating  u between  Eqs.  (19)  and  (20),  the  ratio  of  thrust  to  power 
become  s 


velocity  is  related  to  the  specific  impulse 
= u /g  . (9) 


T 


2 e. 


R 


g 1 


sp 


Note  that  the  time  average  laser  power  is  the  average  continuous  power 
delivered  by  a repetitively  pulsed  device.  The  peak  laser  power  is  the 
average  power  delivered  over  the  duration  of  the  laser  pulse.  The  latter 
concept  is  not  used  in  this  report 
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where  e^  is  defined  as  the  relative  efficiency  and  is  expressed  as 


(10) 


The  quantity  e is  a measure  of  the  energy  efficiency  of  a pulsed  device 

relative  to  a continuous  working  (CW)  or  steady  state  device.  For  CW 

laser  propulsion,  a steady  flow  system  is  obtained  and  the  exit  plane 

velocity  is  the  same  for  all  mass  elements.  Hence,  e is  identically  unity 

and  the  CW  device  is  the  reference  state  for  the  relative  efficiency.  From 

the  mathematical  definition  of  u , it  follows  that 

e 


e 


R 


s 1 . 


To  understand  the  physical  meaning  of  the  relative  efficiency, 

consider  a nozzle  from  which  mass  M/2  exits  with  velocity  V + e and  mass 

M/2  exits  with  velocity  V - S.  The  net  momentum  is  MV  but  the  energy 

2 2 

required  is  1/2  M (V  +e  ).  Hence,  the  thrust  to  power  is  reduced  with 
increasing  non- uniformity  e.  For  pulsed  laser  propulsion,  the  energy 
deposition  is  not  uniform  and  lower  efficiencies  will  result.  The  fluid 
mechanical  model  for  pulsed  laser  propulsion  is  now  extended  to  include 
pulse  sequencing  and  the  relative  efficiency  of  the  pulsed  device  is  assessed. 
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IV.  MUL  ril  I E PI  LSE  THEORY 


I 


The  single  pulse  model  just  treated  assumes  that  the  source  fl  w 
is  established  prior  to  the  detonation  at  r 0.  The  high  pressure  created 
by  the  blast  will  stop  the  source  flow  until  such  time  that  the  pressure  at 
r 0 drops  below  the  pressure  of  the  sonic  orifice.  The  time  at  which  he 
point  source  again  generates  a finite  mass  flux  is  denoted  by  t , which, 
by  definition,  occurs  when 


Ps  (tg  ) g (0)  p 


o r 

/ V 1 /4 


‘s  I */g  (0)  wtf)  r /U  • 

Numerical  results  (presented  later  in  Fig.  6)  indicate  that  t is  signifi- 
cantly shorter  than  the  time  required  for  the  gas  to  flow  through  the 
throat  (r  /u  ).  We  may  w'ish  to  consider  a mechanical  delay  in  re- 
starting the  point  source.  Hence,  we  define  tg  as  the  time  at  which  we 
mechanically  allow  the  source  to  re-start. 


Once  the  point  source  is  re-started,  the  gas  will  expand  with 

velocity  u into  the  relative  vacuum  created  by  the  tail  of  the  blast  wave. 

I 

At  some  time  sufficiently  greater  than  t , we  wish  to  create  a second  pulse. 
Let  us  define  the  critical  time  (tc)  for  the  next  pulse  such  that  the 
shock  would  just  propagate  through  the  source  gas  when  ihe  shock  itself 
becomes  weak.  A schematic  diagram  is  presented  in  Fig.  5a  to  illustrate 
the  meaning  of  tc.  The  value  of  tc  is  readily  determined  by  balancing 


v preceding  paoe 



the  mass  released  by  the  source  during  the  time  interval  (tc  - tg)  with 
the  mass  swept  out  by  the  spherical  shock  as  it  passes  through  the  source 
gas.  Therefore, 


p u (r  ) Q (t  t ) 
c S 


p V n R d t , 
1 s s 


R (t  1 (t  - t ) u 
s max  csx 


Since  t is  defined  by  V (t  ) = u , it  follows  that 

max  s max  X 


tr  - t 


2 max 


We  wish  to  restrict  the  time  between  pulses  (t  ) such  that 

t s t , (12 

pc’  v 

otherwise,  the  shock  will  become  weak  before  it  propagates  through  all  of 
the  source  gas  and  a lower  specific  impulse  will  result.  For  pulse  repeti- 
tion frequencies  which  satisfy  Eq.  (12),  the  shock  will  "break  through" 

the  source  gas  in  a time  (after  detonation)  less  than  t . Let  us  denote 

max 

this  time  by  t^,  as  depicted  in  Fig.  5b. 

We  may  evaluate  t^  in  a manner  equivalent  to  that  used  in 
determining  t . Mass  balance  again  requires 


R (t,  ) = (t  - t ) u , 
s b p s X 
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or 


'b  ('‘n-'.l  U,M\3/2. 


p s £ 


(13b) 


At  time  tfa  (after  detonation)  we  have  a high  pressure  gas  which  must 
expand  isentropically  to  the  exit  plane  of  the  rocket  nozzle.  The  gas  is 
located  between  r = 0 and  r = Rg  (tfa).  The  gas  density,  pressure  and 
velocity  are,  respectively. 


and 


where 


% °sv  n>  • 


pb  p <*b>  «v  ■ 


u.  U (t  ) h(T)  ) 
b s b b 


P|b  = r/RA>  • 


Each  mass  element  dm. 


dm 

“m" 


n 0 (t.  ) R3  (t.  ) f(TL  ) T)  2 dr\ 

s b s b b b b 


- * 2 , 

p u (r  ‘ ) Q (t  - t ) 
P s 


may  expand  to  the  maximum  exit  plane  velocity 


u = u 
e s 


;V  ^ h2(\)  + Yg(Tb)/f(Tb)  , 


corresponding  to  zero  temperature  and  pressure.  The  mass  average  exit 
plane  velocity,  Eq.  (8)  becomes 


2I2  “X 


where 


(y-1 ) J T 


t - t 

JL  s 


(14) 
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and 
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d rj 
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where  we  have  noted  that 
expressions  for  ue  and  u 
tive  efficiency 


the  integrals  over 
, Eq . (10)  may  be 

(Y  + 1)  I/ 


whe  re 


Tl|  and  T)  are  equivalent.  From 
evaluated  to  determine  the  rela- 


(15) 


d ri 


o 

Numerical  results  for  en  along  with  I,  I0,  I and  t u /r'  are  illustrated 

K 1 i.  i s 

in  Fig.  6 and  the  results  are  very  promising.  For  all  practical  purposes, 
the  pulsed  laser  system  is  as  efficient  as  the  CW  system.  The  reason  for 
this  result  is  that  the  bulk  of  the  mass  lies  immediately  behind  the  shock 
where  the  velocity  is  approximately  the  same  for  all  mass  elements.  Hence, 
the  relative  efficiency,  as  defined  by  Eq.  (10),  is  nearly  unity. 

Having  shown  that  the  pulsed  laser  propulsion  concept  is  as  efficient 
as  the  CW  operation,  we  now  use  Eqs.  (14)  and  (2)  to  determine  the  specific 
impulse. 


I 

!£_ 

Vgl 


( y- 1 ) «/ t 


(16) 


The  specific  impulse  is  normalized  to  u ^ / g , the  specific  impulse  of  the 
source  gas  without  energy  addition,  and  the  results  are  illustrated  in  Fig.  7. 
Note  that  t<  1 is  required  in  order  to  obtain  a specific  impulse  significantly 
greater  than  that  of  the  source  gas  alone.  This  is  consistent  with  the 
previously  imposed  constraint  that  t must  be  less  than  t^. 
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Vacuum  Expansion 


Fig.  7 Enhancement  of  Specific  Impulse  Due  to  Pulsed 


Equation  (16)  also  illustrates  that  the  specific  impulse  becomes  unbounded 
as  T -*  0.  This  is  due  to  the  fact  that  the  time  average  laser  power  is  also 
unbounded  in  this  limit.  Rewriting  the  specific  impulse  in  dimensional 
form,  we  obtain 


SP  g 


0 u (D  ) (t  - t ) 
P s 


where  D is  the  throat  diameter  of  the  rocket  nozzle  and  is  related  to  the 
source  radius  r by  continuity, 

* * * 2 * * * 2 

M = 0 u fl  (r  ) = o u TT  (D  ) /4  . 


From  Eq.  (17),  we  conclude 


1.4\/PV/Z 


sp  l g 


where  the  time  average  mass  flow  is  defined  as 

— M (t  - t ) 

M = -M.  = E — , 

t t 

P P 

and  the  time  average  laser  power  is 

P = E/t  . 

P 

The  theoretical  limit  of  the  specific  impulse  is  illustrated  in  Fig.  8.  How- 
ever, these  results  are  general  and  can  be  developed  from  the  general 
relationships  given  in  Sec.  III. 


Vacuum  Expansion 


TIME  AVERAGE  MASS  FLOW 


Subsequently,  Eq.  (17)  is  actually  a general  expression  evaluated  for  that 
value  of  which  corresponds  to  pulsed  laser  propulsion. 

It  is  interesting  to  note  the  scaling  of  the  specific  impulse  with 
molecular  weight.  First,  we  recall  that  in  the  chemical  rocket,  the 

>|c 

specific  impulse  scales  in  direct  proportion  to  u . At  constant  temperature, 
the  specific  impulse  scales  as  the  inverse  square  root  of  the  molecular  weight, 
thereby  favoring  a low  molecular  weight  propellant.  For  the  case  of  pulsed 
laser  propulsion,  we  shall  retain  both  the  chamber  pressure  and  the  laser 

# a{<  sjs  ^ 

power  constant.  Rewriting  Eq.  (17)  with  o replaced  by  y p /(u  ) , we  obtain 


* ' 

u E 

Y p*  <D*>  <tp  - t.l  , 


which  illustrates  that  the  I scales  as  the  square  root  of  u and  the  inverse 

sp 

fourth  root  of  the  molecular  weight,  thereby  favoring  a low  molecular 
weight  propellant,  but  only  by  as  large  a margin  as  in  the  chemical  rocket. 

The  results  illustrated  in  Fig.  8 are  appropriate  only  for  nozzles 
exhausting  with  negligible  temperature  and  pressure.  In  addition,  we  have 
inherently  assumed  that  there  is  no  divergence  of  the  flow  at  the  nozzle 
exit  plane.  The  loss  in  thrust  due  to  a finite  ambient  pressure,  and  that 
due  to  the  finite  size  and  shape  of  rocket  nozzles,  will  now  be  determined  in 


order  to  demonstrate  the  operational  limitations  of  pulsed  laser  propulsion. 


V.  OPERATIONAL.  LIMITATIONS 


I here  are  several  effects  which  limit  the  specific  impulse  that 
we  may  obtain  from  a pulsed  laser  propulsion  system  such  as  that 
analyzed  in  the  previous  sections.  They  are: 

1)  Finite  divergence  of  the  nozzle 

2)  Finite  length  of  the  nozzle. 

3)  Finite  ambient  pressure. 

These  three  effects  can  be  included  in  the  present  fluid  mechanical 

analysis  of  pulsed  laser  propulsion. 

In  determining  the  thrust  from  a pulsed  laser  propelled  rocket, 
we  have  assumed,  without  explicitly  stating,  that  the  rocket  nozzle  is 
contoured  such  that  the  initial  conical  flow  may  exhaust  from  the  nozzle 
without  divergence.  In  the  absence  of  any  such  contouring,  some  thrust 
will  be  lost  due  to  the  divergence  of  the  flow.  The  effect  is  easily 
assessed  by  determining  the  mass  average  of  u^  cos  9,  rather  than  the 
mass  average  of  u . The  thrust  is  reduced  by  the  factor 

2 

sm  ec 

2(1-  cos  o ) 

where  9 is  the  cone  half  angle  and  is  related  to  q by 

n 2 TT  ( 1 - cos  8 ) . 

Results  are  illustrated  in  Fig.  9.  For  cone  half  angles  of  10°,  20°  and 
30°,  the  thrust  is  reduced  by  1%,  3%  and  7%,  respectively.  Hence,  nozzle 


THRUST  (fi>c)  / THRUST  (0C=O) 


8C  - CONE  HALF  ANGLE  (degrees) 

Fig.  9 Thrust  Reduction  Due  to  Fiow  Divergence 
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intouring  will  be  a refinement,  rather  than  a primary  fluid  mechanical 
consideration.  A preliminary  assessment  of  the  effects  of  a continuously 
changing  nozzle  wall  angle  (eg.  a parabolic  nozzle)  is  presented  in  the 
Appendix. 

The  length  of  the  nozzle  is  restricted  by  the  radius  of  the  shock  at 
'breakthrough1'.  The  shock  must  be  contained  within  the  nozzle,  other- 
wise the  specific  impulse  will  be  reduced.  This  is  expressed  as 


L > R (t  ) 
s b 


(18) 


where  L is  the  nozzle  length  and  R (t^  ) is  the  shock  "break  through1 
radius,  given  by  Eq.  (13). 


R (t.  ) 

s b 


t - t 

_E ,s’ 

d"7u" 


D 


(19) 


However,  satisfying  Eq.  (18)  does  not  insure  that  the  maximum  thrust  is 

obtained  from  the  energy  deposition.  In  order  to  obtain  the  results  quoted 

in  Fig.  8,  the  gas  must  expand  to  a vacuum  while  still  within  the  nozzle. 

To  determine  the  effect  of  the  finite  value  of  L on  the  I , we  allow  the 

sp 

gas  to  expand  i sentropically  to  the  exit  plane  at  r - L.  We  conserve  mass, 
entropy  and  energy  by 


2 2 2 
D u L - o.  u r\  R (t  ) 
e e b b b s b 


(20) 


and 


P,  / 0 ^ p /o  Y 
b b e e 


(21) 


C T + 1/2  u = C T,  + 1/2 
p e e p b 


(22) 


-31- 


L 





respectively,  where  the  subscripts  b and  e denote  "break  through"  and  exit 


plane  conditions,  respectively.  Solving  for  u , we  obtain 

e 


ue  = Us(tb)  V H (V  + Y YCy  8 (\)/f  (\>  . (23) 


where  Y (T^)  is  a weak  function  of  u and  must  be  determined  by  iteration. 


Denoting  the  Nth  iteration  by  Y 1 , we  obtain 


y<n>  = i . (T1  /*)2<Y-n  ( u* 

1 V ' l (N-l) 


where 


1 = L/Rs(tb)  , 


u = h on  ) U (t  ) , 

e b s b 


(N)  (N) 

and  u^  is  determined  from  Y via  Eq.  (23)  for  N 2=  1. 


The  reduced  thrust  and  specific  impulse  due  to  the  finite  nozzle 


length  is  obtained  by  using  the  corrected  value  of  u in  the  integrand  of  I . 

e 2 


Results  are  illustrated  in  Fig.  10.  For  y = 1.4,  termination  of  the  nozzle 
at  the  break  through  radius  would  result  in  a recovery  of  only  68%  of  the 
thrust.  However,  90%  of  the  thrust  is  recovered  in  three  break  through 


radii.  Nozzle  lengths  greater  than  three  (t^)  yield  a deminishing  return 


in  terms  of  thrust  recovered.  We  shall  use  L = 3 R (t,  ) as  a design 

s b 


criterion. 


The  finite  ambient  pressure  will  also  reduce  the  specific  impulse 
of  the  laser  powered  rocket.  To  illustrate  the  magnitude  of  this  effect, 


.J 


we  again  allow  the  gas  to  expand  isentropically  from  the  "break  through" 
conditions  to  the  ambient  pressure,  p^.  Conservation  of  energy  and 
entropy,  Eqs.  (22)  and  (21),  yield 
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and 
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where  Q is  approximately  the  ratio  of  p^  to  the  pressure  of  the  blast  wave. 

The  reduced  thrust  and  specific  impulse  due  to  the  finite  ambient 

pressure  is  obtained  by  using  the  corrected  value  of  u^  in  the  integrand  of 

I . Results  are  illustrated  in  Fig.  11.  For  a chamber  pressure  of  two 
2 

atmospheres  (sufficient  to  "choke"  the  flow  at  sea  level),  the  sonic  pressure 
is  one  atm.  Laser  energies  sufficient  to  create  a 30,  000  I<  plasma  would  result 
in  a blast  wave  pressure  of  100  atm.  Thus,  the  sea  level  value  of  would  be 
10”^  and,  for  y = 1.4,  75%  of  the  vacuum  thrust  would  be  obtained.  Raising  the 
chamber  pressure  to  20  atm  would  result  in  a 1000  atmosphere  blast  wave  and 
a recovery  of  90%  of  the  vacuum  thrust.  Higher  chamber  pressures  yield 
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diminishing  returns  in  the  thrust  recovered  and  would  begin  to  introduce 
thermodynamic  and  structural  penalties. 

Having  determined  the  limitations  for  which  the  results  of  Sec.  IV 
are  valid,  we  will  illustrate  the  results  by  way  of  a simple  example.  From 
Sec.  Ill,  we  recall  the  relation  between  power  P and  thrust  T 

P = g I TV  2 , (25) 

sp 


where 


R 


1 and 


P = E/t  . (26) 

P 

Equations  (17)  and  (25)  yield  the  size  of  a rocket  required  for  the  thrust  T. 


If  we  wish  to  design  a rocket  to  operate  with  a chamber  pressure  of  two 

atmospheres  and  temperature  of  300°K,  the  sonic  density  and  velocity 

become  10"3  gm/cc  and  30,  000  cm/sec,  respectively.  With  a specific 

impulse  of  800  seconds,  we  may  consider  launching  a one  ton  rocket  with 

10  g' s acceleration.  This  requires  an  orifice  diameter  of  22  cm  and  a 

power  of  350  megawatts,  as  indicated  in  Figs.  12  and  13,  respectively. 

In  a laboratory  experiment  (to  be  discussed  in  Sec.  VII)  we  may  consider 

an  equivalent  acceleration  of  a 2 pound  rocket  which  requires  D = 7 mm  and 

P .35  megawatts,  also  indicated  in  Figs.  12  and  13  respectively.  The 

remaining  point  to  emphasize  is  that  there  are  distinct  combinations  of  E and 

t permitted  by  the  nozzle  fluid  mechanics.  These  are  illustrated  in  I ig.  14. 
P 
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where  we  illustrate,  from  Eq.  (17),  t vs  D at  constant  E,  For  arbitrary 

P 

values  of  the  orifice  diameter,  D , there  is  an  upper  limit  on  the  pulse 
repetition  frequency  for  which  the  point  source  approximation  is  valid. 

This  reflects  the  fact  that  a finite  time  is  required  for  the  source  gas  to 
enter  the  rocket  nozzle.  There  also  exists  a lower  limit  on  the  pulse 
repetition  frequency  for  which  a nozzle  of  a given  length  can  recover  the 
optimum  thrust  from  the  energy  deposited.  Larger  time  durations  between 
pulses  would  permit  the  source  gas  to  escape  the  nozzle  before  all  of  the 
pressure  were  converted  to  momentum.  Our  design  criterion  required 
that  the  nozzle  length  be  three  "break  through"  radii  to  recover  90%  of  the 
thrust.  Nozzle  lengths  greater  than  50  D are  excluded  on  the  basis  of 
their  large  aspect  ratio.  The  illustrative  examples  indicate  that  the 
laboratory  experiment  could  be  conducted  with  a 15  joule  laser  operating  at 
25,  000  pulses  per  second.  Launching  one  ton  with  20,000  lbf  of  thrust 
would  require  a megajoule  laser  operating  at  350  pulses  per  second.  The 
conditions  for  both  the  laboratory  experiment  and  the  "full  scale"  rocket 
are  well  within  the  realm  of  possibility. 

Finally,  using  Figs.  12  - 14,  we  can  determine  the  combination  of 
laser  energy  and  pulse  repetition  frequency  (and  therefore,  average  laser 
power)  required  to  obtain  a specific  impulse  of  800  sec.  with  a pulsed  laser 
propulsion  system.  In  Fig.  15  laser  energy  versus  pulse  repetition 
frequency  (PRF),  along  with  lines  of  constant  average  power  (equal  to 
E/tp),  is  presented.  The  I = 800  sec.  design  curve,  obtained  from  Figs  . 

12  - 14,  is  shown.  In  addition,  a recommended  laser  pulse  time  for  the 
specified  laser  energy  is  presented  on  the  right  hand  ordinate.  This  scale 
is  obtained  using  blast-wave  theory  along  with  the  requirement  that  the 
blast  wave  does  not  propagate  greater  than  one -half  the  throat  diameter 
during  the  duration  of  the  laser  pulse. 
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Fig.  15  Laser  Energy  vs.  Pulse  Repetition  Frequency 
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LASER  PULSE  TIME  (sec) 


VI.  NOZZLE  AND  PLENUM  PR ESSUR  E TRACES 


Practical  engineering  considerations  make  it  necessary  to  illustrate 
the  pressure  distribution  in  the  nozzle  and  plenum  chamber  (from  which 
the  propellant  is  supplied)  as  a function  of  space  and  time.  The  pressure 
behind  the  blast  wave  is  given  in  Sec.  II. 


ps  (t)  (ttt) 

Using  Eqs . (1)  and  (4)  we  obtain 
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which  is  illustrated  in  Fig.  16.  For  E - one  megajoule,  y - 1.4  and 

3 

1,  the  peak  pressure  is  about  10  atm  at  the  throat,  decreasing  to 
below  one  atm  at  L = 3 m. 

The  time  history  of  the  pressure  at  the  throat  is  determined  from 


Pg  (t)  and  g (0),  obtained  from  Fig.  3. 


P (0,  t) 


1/2  .,.2 

P (D ' ) g (0) 


This  pressure  time  history  is  illustrated  in  Fig.  17  for  the  one  ton  rocket. 
The  throat  pressure  relaxes  to  below  p in  10"^  seconds,  which  is  small 

5 » 

with  respect  to  D /u  . Therefore,  the  ''source"  restarts  and  establishes  a 
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Pressure  Time  History  in  Throat 
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new  nozzle  flow  before  the  next  pulse  initiates  another  blast  wave.  The 
results  quoted  in  Fig.  17  should  be  treated  with  caution.  They  are  the 
results  for  an  instantaneous  point  source  of  energy.  We  are  using  these 
results  in  spatial  locations  that  are  of  the  order  of  the  deposition  zone,  and 
therefore  some  error  is  expected.  Furthermore,  and  most  important,  we 
are  quoting  times  that  may  be  small  with  respect  to  the  energy  deposition 
time.  If  the  laser  pulse  duration  is  greater  than  2 microseconds,  the  early 
time  portion  of  the  pressure  pulse  is  in  error.  At  times  large  with  respect 
to  the  pulse  duration  time,  the  point  source  approximation  becomes  valid. 

A further  complication  in  the  pressure  time  history  at  the  throat  is 
due  to  the  fact  that  the  throat  diameter  is  finite.  Hence,  the  blast  not  only 
propagates  down  the  nozzle,  but  into  the  plenum  chamber  as  well.  This  high 
pressure  plenum  gas  will  also  create  deviations  from  the  predictions  made 
in  Fig.  17.  Suppose  we  model  the  plenum  shock  as  a spherical  blast  wave 
propagating  into  a uniform  atmosphere.  Conventional  blast  wave  theory 
yi  elds 


and 


where  = CQ  (y)  **  1 • The  pressure  at  the  throat  is  approximately  one  half 

of  PSj  and  this  approximation  for  the  pressure  time  history  at  the  throat  is  also 

indicated  in  Fig.  17.  The  plenum  model  yields  a much  higher  pressure  which 

-4 

indicates  that  the  mass  flux  begins  earlier  than  10  seconds  but  the  pressure 
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does  not  relax  to  p until  10"3  seconds.  This  is,  however,  still  well 
in  advance  of  the  next  pulse.  The  location  of  the  plenum  shock  at  the 
time  of  pressure  relaxation  is  also  obtained  from  the  spherical  blast  wave 
model  and  corresponds  to  one  meter  at  10'3  seconds.  This  length  scale 
is  in  excess  of  the  throat  diameter  (22  cm)  and  may  be  as  large  as  the 
plenum  chamber  itself.  This  creates  two  potential  problems.  First: 
each  pulse  will  be  detonated  in  a gas  which  has  already  been  processed 
by  the  plenum  shock  and  secondly,  the  mechanical  pumps  supplying  the 
expellant  to  the  plenum  chamber  may  have  to  work  against  a higher  back 
pressure.  An  assessment  of  these  effects  and  any  correction  thereof  is 
beyond  the  scope  of  this  study. 
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VII.  EXPERIMENT  DESIGN 


An  experiment  designed  to  test  the  pulsed  laser  propulsion  concept 
cun  be  performed  using  pulsed  CO,,  TEA  lasers.  The  objectives  of  such 
an  experiment  are  to  (1)  verify  the  theoretical  predictions  of  high  specific 
impulse  with  relative  efficiency  near  unity,  (2)  examine  the  dependence 
of  the  time  averaged  specific  impulse  and  thrust  upon  average  laser  power 
and  repetition  rate,  (3)  study  the  effects  of  finite  exit  plane  pressure  upon 
thrust  and  specific  impulse  and  (4)  determine  how  the  effects  of  energy  lost 
to  the  walls  along  with  the  finite  source  nature  of  the  breakdown  process 
will  affect  rocket  performance  and  comparison  with  theory. 

Multiple  independently  triggered  pulsed  CC>2  TEA  lasers  can  be 
utilized  for  such  experiments.  The  pulse  shape  of  such  a laser  is  typically 
composed  of  a 100  nsec  spike  with  a one  to  three  gtsec  tail.  The  tail  of  the 
pulse  can  be  virtually  eliminated  by  reducing  the  amount  of  nitrogen  in  the 
active  medium.  Therefore,  the  effects  of  the  tail  of  the  pulse  on  the  rocket 
performance  can  be  examined.  A "fifteen  joule"  CO£  TEA  laser  will 
deliver  5 J within  a 100  nsec  laser  pulse,  or  1 5 J within  a 1 00  nsec  spike 
and  a 3 psec  tail.  The  minimum  number  of  lasers  required  to  achieve 
complete  simulation  of  multiple  pulse  pertormance  is  not  clear.  However, 
at  least  two  pulses  are  necessary  to  determine  if  the  independence  of  each 
pulse,  assumed  theoretically,  is  realistic.  A schematic  of  an  experiment, 
which  utilizes  four  lasers  of  1 5 J each,  is  presented  in  Fig.  18.  The  four 
laser  beams  will  be  directed  into  the  vacuum  chamber  through  a window  at 
slightly  different  angles  of  incidence  (up  to  10-2  radians  apart)  so  the 
focused  spots  will  be  on  the  order  of  a centimeter  in  diameter.  Therefore, 
the  performance  should  be  equivalent  to  that  obtained  with  co-linear  laser 
beams.  The  lasers  are  individually  triggered,  and  thus,  various  interpulse 
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times  can  be  programmed  to  obtain  the  desired  pulse  repetition  rates. 

At  a repetition  rate  of  10,  000  pps,  this  laser  system  will  deliver  an 
average  laser  power  of  150  kilowatts.  The  si/.e  of  'he  rocket  nozzle  and 
the  anticipated  performance  of  such  a system  will  be  determined  from 
theoretical  predictions,  and  will  be  discussed  later. 

Diagnostics 

The  thrust  and  specific  impulse  are  the  two  parameters  which  will 
characterize  the  performance  of  the  pulsed  laser  propulsion  system.  The 
relationship  between  thrust,  average  power  and  specific  impulse  is  given 
by  Eq.  (25),  and  the  average  power  as  a function  of  laser  energy  and 
repetition  rate  is  given  by  Eq.  (26).  The  thrust  anticipated  from  the  four 
15  .1  laser  pulse  system  as  a function  of  specific  impulse  and  repetition  rate 
is  presented  in  Fig.  19.  The  thrust  can  be  measured  either  by  using  a load 
cell  to  obtain  the  instantaneous  thrust  or  by  using  a ballistic  pendulum  to 
measure  the  delivered  total  impulse/pulse  and  dividing  it  by  the  time 
between  pulses  to  determine  the  time  averaged  thrust.  The  ballistic 
pendulum  is  a simple  device  which  directly  measures  the  total  delivered 
impulse.  However,  since  the  pendulum  must  be  free  to  swing,  the  propellant 
must  be  carried  as  an  integral  part  of  the  pendulum  structure.  The  load 
cell  provides  instantaneous  thrust,  but  the  inherent  response  time  of  the 
cell  and  its  response  to  acoustic  ringing  in  the  nozzle  wall  may  make 
qualitative  results  difficult.  It  is  most  desirable  to  use  the  load  cell 
combined  with  nozzle  wall  pressure  transducers  (as  in  Ref.  4)  to  deduce 
the  pres su re /time  history  along  with  the  thrust/time  history.  Comparison 
with  integrated  thrust  (or  impulse)  from  a ballistic  pendulum  will  serve  as 
a calibration  of  the  load  cell  response. 

The  specific  impulse  is  difficult  to  measure  directly.  It  can  be 
obtained  if  the  impulse  and  the  amount  of  gas  expelled  from  the  nozzles  is 
measurable.  A more  promising  approach  to  determining  specific  impulse 
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is  to  directly  measure  the  exhaust  velocity.  The  technique  that  will  be 
used  involves  multiple  sparks  in  the  flow  field  from  the  same  set  of 
electrodes.  * ^ The  hot,  partially  ionized  gas  flows  with  the  instantaneous 
p exhaust  velocity.  Therefore,  after  one  spark  is  ignited,  the  hot  gas 

propagates  at  the  exhaust  velocity.  A second  spark  is  then  ignited  at  the 
new  position  of  the  partially  ionized  gas.  By  measuring  the  distance  between 
the  first  and  second  spark  positions,  the  exhaust  velocity  is  determined 
(see  Fig.  20  and  Ref.  15).  For  CW  measurements  with  a streak  or  framing 
camera,  a version  of  Jacob's  ladder  could  be  made.  The  familiar  Jacob's 
ladder  spark  rises  because  of  free  convection,  while  in  the  present  case 
the  spark  convects  downstream  at  the  exhaust  velocity.  A streak  camera 
measurement  yields  information  along  a specified  section  of  the  flow  field 
while  a framing  camera  is  used  to  obtain  velocity  as  a function  of  distance 
from  the  nozzle  centerline.  However,  neither  of  these  techniques  will 
result  in  the  velocity  at  the  exhaust  plane  as  a function  of  time. 

Nozzle  Design  and  Propellant  Supply  Considerations 

It  is  desirable  to  perform  experiments  in  pulsed  laser  propulsion 
using  two  distinct  nozzle  configurations.  Since  theoretical  modeling  has 
been  completed  for  a conical  nozzle,  it  is  essential  to  perform  a series  of 
experiments  with  such  a nozzle  to  verify  the  theoretical  predictions.  This 
configuration  is  gasdynamically  the  simplest,  and  beam  focusing  is  accom- 
plished externally  which  permits  control  and  variation  of  spot  size  and 
breakdown  location.  On  the  other  hand  since  an  actual  system  may  utilize 
a self-focusing  nozzle,  experiments  with  a parabolic  nozzle  should  also 
be  performed.  This  nozzle  is  optically  optimum  and  permits  near 
spherical  focusing  of  a collimated  laser  beam.  Experiments  with  this 
nozzle  shape  are  necessary  to  assess  any  aerodynamic  penalties  resulting 
from  wall  curvature  (due  to  the  unsteady  nature  of  (he  flow  as  discussed  in 
the  Appendix)  and  determine  the  effectiveness  of  ’he  self  focusing  concept 
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Fig.  20  Example  of  Spark  Technique  for  Measuring  Exhaust 
Velocity  from  Ref.  15 
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to  obtain  a strong  laser  induced  breakdown  at  the  focal  point.  Nozzle 
sizing  is  accomplished  by  using  the  theoretical  predictions  for  optimum 
performance.  Essential  experiments  to  be  performed  are  as  follows: 

1.  Using  a conical  nozzle  in  a vacuum  chamber  with  an 

external  propellant  source,  perform  a series  of  experiments 

to  determine  the  specific  impulse  as  a function  of  average 

laser  power/mass  flow  (Fig.  8).  In  addition,  it  is  desirable 

to  explore  the  operating  corridor  in  repetition  rate  at  constant 

specific  impulse  (e.g.  Fig.  14,  where  I 800  sec). 

s p 

2.  Using  the  conical  nozzle  at  a specific  design  point, 
examine  the  effect  of  finite  exhaust  plane  pressure  on  the 
observed  thrust  and  specific  impulse. 

3.  The  above  experiments  with  two,  three  and  four  laser 
pulses  should  be  performed  to  determine  the  independence  of 
each  thrusting  sequence  along  with  the  pre s s u re /time  history 
in  the  propellant  plenum  chamber. 

14.  A selected  number  of  the  above  experiments  with  a 

parabolic  nozzle  should  be  performed  to  determine  the  focusing 
ability  of  the  nozzle  and  the  aerodynamic  penalties  of  the  wall 
curvature. 

When  sizing  the  nozzle  for  vacuum  experiments,  we  want  to  be 
sure  the  nozzle  is  sufficiently  long  for  the  experiments  which  are  to  be 
performed  at  various  repetition  rates.  Theoretically,  we  have  determined 
that  it  is  desirable  to  have  the  nozzle  length  equal  to  at  least  3Rb  where 
is  the  "breakthrough"  radius.  Equation  (19)  relates  this  radius  to  the 
time  between  laser  pulses  and  the  sonic  velocity  of  the  propellant  at  the 
entrance  to  the  nozzle.  If  the  plenum  is  at  room  temperature  and 
y 1.4,  u « 3. 22x10^  cm  / sec.,  and  from  Eq.  (19) 

R (t  ) ~ 8 x 10  (tp  - tg). 

bp  r 


\ 
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For  lasers  of  5 to  1 5 J it  is  desirable  to  operate  in  a repetition  rate  regime 


such  that  1 0*5  < t - tg  < 10“  . Since  in  a vacuum  there  is  no  penalty  for 

having  the  nozzle  too  long,  it  is  best  to  design  around  (t  - t ) = 10“^. 

P s 


Therefore  (t  ) ~ 8 cm,  and  a recommended  nozzle  length  would  be 
L 24  cm.  Since  it  is  desirable  to  "choke'1  the  propellant  flow  at  sea 


level  for  earth-based  missions  we  shall  calculate  the  plenum  chamber 
conditions  for  an  experiment  based  upon  p 1 atm.  The  resulting  stagna- 


tion pressure  for  y 1 . 4 becomes  pQ  = 1.9  atm  and  with  a plenum 
temperature  TQ  300°K,  we  obtain  p 2.  13  x 10  ^ gms/cm^  and 
p 1.35  x 10  gms/cm  . In  order  to  determine  the  propellant  mass  flow 


rate  it  is  necessary  to  choose  a throat  diameter.  The  best  way  to  operate 
in  the  corridor  of  Fig.  14  is  fix  the  throat  diameter  and  vary  the  laser 


energy  and  the  time  between  pulses.  For  laser  energies  between  5 and  15  J 


a throat  diameter  D'1’  = 0.  5 cm.  will  permit  operation  at  I = 800  sec. 

5 S ^ 4 

when  the  time  between  laser  pulses  is  2 x 10"  < tn  - t < 10  sec. 

2 S 

If  D = 0.  5 cm.  , A ^ . 2 cm  and  the  mass  flow  rate  equals  8.  5 gms  /sec. 


The  volume  of  the  propellant  supply  chamber  (or  plenum)  required  to 
perform  an  experiment  with  four  laser  pulses,  for  example,  can  now  be 


calculated.  If  the  lasers  are  pulsed  at  the  slowest  repetition  rate  of 


1 04  pps  and  the  la'ser  pulse  time  is  2 usee,  the  entire  experiment  is 


completed  in  approximately  4 x 10  sec.  At  a propellant  mass  flow  of 


8.5  gms/sec.  this  corresponds  to  a total  propellant  mass  utilized 


~ 3 -3  3 

M ~ 3.5x10“  gms . If  pQ  = 2,  13  x 10  gms  / cm  , the  volume  of 


propellant  utilized  fiV  = 1.6  cm  . It  is  desirable  to  have  the  change  in 


propellant  volume  not  exceed  10%  of  the  total  propellant  volume.  There- 


fore, a plenum  of  16  cm  is  desired. 


A schematic  of  the  above  nozzle  designed  for  experiments  with  four 
pulsed  TEA  lasers  of  5 to  15  J is  shown  in  Fig.  21.  Since  the  unit  should 
be  completely  self  contained  for  operation  in  a vacuum  chamber,  the 


n 
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Fig.  21  Conical  Nozzle  Design 


plenum  is  attached  to  (he  nozzle  with  a diaphragm  inserted  to  maintain 
the  stagnation  pressure  in  Ihe  plenum  until  the  experiment  is  performed. 
Precise  timing  of  the  bursting  of  the  diaphragm  prior  to  the  laser  pulse 
is  desirable  in  order  to  prevent  the  vacuum  chamber  from  filling  up  with 
propellant,  raising  the  chamber  pressure.  With  a supersonic  propellant 
flow  in  the  nozzle  the  exhaust  plane  pressure  should  be  maintained  at 
2 x 10  ^ torr  since  the  ratio  of  the  exit  plane  area  to  the  throat  area  is 
approximately  700.  Use  of  a diaphragmed  plenum  may  make  turn-around 
time  between  experiments  longer  than  desired,  but  experiments  with  a 
steady  source  of  propellant  into  the  nozzle  would  require  a large  pumping 
capacity  if  experiments  are  to  be  performed  in  a vacuum. 

In  Fig.  2 1 the  nozzle  half  angle  has  been  taken  to  be  15°  . Thus, 
from  Fig.  9 greater  than  95%  of  the  theoretical  thrust  estimate  can  be 
anticipated  with  this  nozzle.  This  angle  also  does  not  impose  severe 
restrictions  upon  the  beam  focusing  requirements  which  are  discussed 
in  the  next  subsection.  Also  shown  in  Fig.  21  are  the  suggested  optimum 
locations  for  wall  mounted  pressure  transducers.  Transducer  No.  1 is 
mounted  in  the  plenum  wall  in  order  to  monitor  the  change  in  the  stagnation 
pressure  over  the  four  thrusting  pulses.  Transducer  No.  2 is  mounted 
as  close  to  the  throat  as  possible  in  order  to  obtain  a measurement  of  the 
pressure  in  the  breakdown  region.  Transducer  No.  3 is  placed  at  the 
estimated  "breakthrough"  radius  to  determine  the  propellant  pressure 
prior  to  the  isentropic  expansion,  and  transducer  No.  4 is  mounted  in  the 
exhaust  plane. 

A schematic  of  a parabolic  nozzle  arrangement  is  shown  in  Fig.  22 
This  nozzle  was  designed  such  that  A^/A  , where  A^  is  the  exit  plane 
area,  is  the  same  as  the  cone  nozzle,  and  the  angle  at  which  the  flow 
leaves  the  nozzle  is  15°,  as  it  is  with  the  cone  nozzle.  The  corresponding 
length,  however,  is  now  only  1.5  "breakthrough"  radii.  The  overall  effect 
of  the  shortened  length  and  parabolic  shape  of  the  nozzle  must  be  assessed 
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experimentally.  Again,  suggested  pressure  transducer  locations  are 
indicated  in  Fig.  22  at  the  plenum  wall,  the  throat,  the  "breakthrough" 
radius  and  the  exit  plane. 

Laser  Beam  Focusing  Requirements 


The  most  important  focusing  requirement  is  that  a laser -induced 

breakdown  is  achieved  in  the  vicinity  of  the  throat.  For  most  of  the 

, 9 . 2 

anticipated  propellants  a focused  laser  intensity  of  > 1 0 W/cm  will  yield 

a breakdown  in  the  1 atm  gas  at  the  throat.  The  peak  power  from  the 

above  mentioned  lasers  can  be  computed  from  the  amount  of  laser  energy 

in  the  spike.  Using  5 J in  a 100  nsec  pulse,  we  obtain  a peak  power  of  50  MW. 

For  the  conical  nozzle  in  which  the  beam  is  focused  externally,  mirrors 

-3  2 

can  focus  each  individual  beam  to  a spot  area  of  10  cm  . Therefore, 
the  peak  intensity  at  the  focus  will  exceed  10^  W/cm^.  Achieving  a break- 
down should  not  be  a problem  for  conical  nozzle  experiments.  The  focusing 
ability  of  a parabolic  nozzle  will  depend  upon  many  factors  including  the 
surface  condition  and  the  precise  curvature  of  the  nozzle  walls.  The 
external  optics  becomes  simplified,  and  a series  of  flats  can  be  used  to 
pass  each  beam  into  a nozzle  (see  Fig.  18).  The  strength  of  the  laser- 
induced  breakdown  is  one  of  the  desired  measurements  to  be  made  using  a 
wall  mounted  pressure  transducer. 

Propellant 

The  choice  of  a propellant  will  be  dictated  by  (1)  the  desire  to 
obtain  a high  specific  impulse,  (2)  ease  of  obtaining  a breakdown  and 
(3)  convenience  for  laboratory  operation.  Since  it  is  desirable  to  use  a 
low  molecular  weight  propellant,  most  experiments  can  be  performed 
using  helium.  A few  experiments  with  hydrogen  should  be  attempted 
since  ultimately  hydrogen  may  be  more  desirable  from  a systems  view- 
point. The  use  of  other  propellants  (such  as  water  vapor)  that  may  prove 
advantageous  as  a result  of  systems  analyses  should  be  reviewed. 
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Test  Plan 


A typical  test  plan  for  a possible  pulsed  laser  propulsion  experi- 
ment is  presented  in  Fig.  23.  Before  embarking  on  a series  of  experiments 
with  multiple  laser  pulses  and  propellant  flow,  however,  it  is  desirable 
to  perform  a series  of  single  pulse  experiments  with  only  ambient  air  in 
the  nozzle,  similar  to  the  experiments  performed  in  Ref.  4.  It  was  found 
in  Ref.  4 that  long  laser  pulses  lead  to  inefficient  conversion  of  laser  power 
to  thrust  with  a resulting  low  specific  impulse.  Since  the  TEA  lasers  have 
short  pulses,  a series  of  experiments  with  no  propellant  flow  should  be  per- 
formed to  determine  the  dependence  of  the  propulsion  efficiency  on  pulse 

time.  These  experiments  should  be  performed  at  various  laser  energy 

1 3 

levels  and  exit  plane  pressures.  Back  has  presented  a theory  for  the 

specific  impulse  expected  from  a single  blast  at  the  apex  of  a conical  nozzle. 

Comparisons  with  this  theory  for  the  conical  nozzle  experiments  can  be 

made.  In  addition,  such  experiments  can  be  used  to  calibrate  a load  cell 

and  examine  the  focusing  ability  of  the  parabolic  nozzle  configuration. 

The  test  matrix  in  Fig.  23  has  been  divided  into  two  parts.  The 

first  series  of  experiments,  where  specific  impulse  is  the  output,  is 

designed  to  examine  the  validity  of  Fig.  8 and  Eq.  17.  For  the  experimental 

conditions  enumerated  in  the  previous  section  each  combination  of  laser 

energy  (E)  and  time  between  pulses  (t  - tc)  will  yield  a time  averaged 

delivered  power,  and  using  Eq.  (17)  or  Fig.  8,  we  calculate  the  expected 

specific  impulse  (I  ) presented  in  the  matrix.  Experiments  should  be 
sp 

performed  for  this  variety  of  laser  conditions  first  in  a vacuum,  then 

with  increasing  exit  plane  press  .re,  and  with  both  the  conical  nozzle  and 

the  parabolic  nozzle.  The  second  series  of  experiments,  where  we  want  to 

examine  the  constant  I corridor  of  Fig.  14,  is  designed  to  determine  the 

sp 

regime  of  validity  of  the  theoretical  model.  These  experiments  should  be 
performed  with  the  laser  energy  and  time  between  pulses  that  is  specified 
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experiments  to  be  performed  with  conical  and  parabolic  nozzles. 


( SEC  ) 


to  yield  a constant  specific  impulse.  Again,  the  experiments  should  be 
performed  with  both  nozzles  and  at  various  exit  plane  pressures  to 
determine  the  vacuum  and  atmosphere  performance  to  be  expected  from 
the  pulsed  laser  propulsion  system. 


<! 

I 
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VIII.  SUMMARY  AND  CONCLUSIONS 


A fluid  mechanical  model  to  describe  the  flow  within  a conical  nozzle 

that  is  subjected  to  point  energy  depositions  at  the  apex  of  the  cone  has  been 

developed.  The  model  has  been  used  to  assess  the  concept  of  pulsed  laser 
(4) 

propulsion.  The  specific  impulse  of  a pulsed  laser  propelled  rocket  has 
been  obtained  as  a function  of  nozzle  sonic  conditions,  the  laser  energy  and 
the  pulse  repetition  frequency.  The  results  are  given  by  Eq.  (17)  and 
illustrated  in  Fig.  8.  Specific  impulses  of  the  order  of  a few  thousand 
seconds  may  be  obtained  with  moderate  energy  densities  and  pulse  repeti- 
tion frequencies.  However,  these  results  are  valid  only  for  an  expansion  ] 

to  a vacuum.  The  reduction  in  the  thrust  due  to  a finite  exit  plane  pressure 
may  be  significant.  The  finite  exit  plane  pressure  may  arise  due  to  either 
the  finite  ambient  pressure,  or  the  finite  length  of  the  rocket  nozzle.  The 
degradation  of  the  specific  impulse  due  to  both  of  these  effects  has  been 
assessed. 

The  relative  efficiency  (as  defined  in  Eq.  10)  of  a pulsed  propulsion 
system  has  been  determined  and  compared  to  a C W system.  This  efficiency 
is  based  on  fluid  mechanical  considerations  only,  and  the  laser  absorption 
mechanisms  have  not  been  considered.  That  is,  we  have  assumed,  for 
both  CW  and  pulsed,  that  all  of  the  laser  energy  has  been  absorbed  by  the 
working  medium.  Results  indicate  that  pulsed  laser  propulsion  is 
approximately  98%  to  99%  as  efficient  as  the  CW  device  in  converting 
laser  power  to  thrust.  Hence,  pulsed  laser  propulsion  appears  as  versatile 
and  as  efficient  as  the  CW  technique  with  the  potential  advantages  of 
simplicity  in  engine  design  and  the  elimination  of  possible  plasma  stability 
constraints  that  may  be  associated  with  a CW  laser  propulsion  system.  Finally, 
an  experiment  to  test  the  theoretical  model  and  validate  the  overall  concept 
has  been  designed  and  appears  feasible  using  existing  commercially  produced 
la  s e r s . 
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APPENDIX 


SHOCK  REFLECTION  FROM  NOZZLE  WALLS 


In  analyzing  the  pulsed  laser  propulsion  concept,  we  have  assumed 
that  the  nozzles  are  conical.  In  reality,  the  nozzles  will  have  a shape  that 
is  other  than  conical  and  the  blast  wave  will  reflect  from  the  nozzle  walls. 
These  shock  reflections  will  create  pressure  and  velocity  perturbations  in 
the  flow  which  will  degradate  the  performance  of  the  system.  This  degrada- 
tion must  be  assessed  in  order  to  insure  that  the  pulsed  system  is  fluid 
mechanically  as  efficient  as  the  CW  device. 

For  simplicity,  we  will  consider  the  "real  nozzle"  illustrated  in 

Fig.  A- 1.  The  nozzle  is  assumed  to  be  conical  for  length  R , where  the 

o 

curvature  is  concentrated  into  a single  turn  in  the  wall.  The  deflection  of 
the  wall  creates  a corner  shock,  A,  to  turn  the  flow  parallel  to  the  wall,  and 
a lamda  shock,  B,  which  is  the  reflection  of  the  blast  wave  from  the  wall. 
The  shock  structure  will  propagate  radially  toward  the  nozzle  centerline  as 
the  blast  wave  moves  down  the  nozzle.  Ultimately,  the  shock  structure  will 
reflect  from  the  nozzle  centerline  and  form  a second  normal  shock,  further 
complicating  the  flowfield.  Our  problem  is  to  first  determine  that  portion 
of  the  flowfield  which  will  cause  the  primary  thrust  degradation. 

Figure  A-2  illustrates  that  portion  of  the  thrust  which  is  generated 
by  the  gas  between  T|  0,  and  TU  Clearly,  90%  of  the  thrust  is  generated  by 
the  gas  between  1 0,6  and  )=  1.  Hence,  we  are  concerned  only  with  shock 


A-  1 


THUST  BETWEEN  0 AND 


reflections  which  occupy  that  portion  of  the  gas.  We  will  demonstrate  that 
the  second  normal  shock  cannot  form  in  the  T|  > 0.  6 region,  and  the  primary 
fluid  nechanical  requirement  is  to  model  the  flowfield  illustrated  in  Fig.  A-l. 

Region  1 of  Fig.  A-l  is  a zone  of  silence.  The  fluid  mechanical 
properties  of  region  1 have  been  obtained  in  Sec.  III.  The  shock  radius 

R (t),  increases  as 

s 

2/3 

R (t)  = At  1 , 

s 

where  A is  a constant  of  the  motion.  The  fluid  velocity  is 


u 


(ft)  vs  h m ■ 


where 

V = (2/3)  A t 
s 

and 


h (T|)  « T1  , 


where 


•H  = r/Rg  (t) 


The  Mach  number  of  the  gas  is  denoted  by 


M j - u / yRT 


(A-l) 


where  we  have  made  the  approximation  that  the  temperature  is  constant  in 
the  shock  profile.  Numerical  results  indicate  that  Eq.  (A-l)  is  within  5% 
accuracy  for  y = 1.  4 and  0.  2 < *1  < 1. 


A-4 


i'li e relatively  low  Mach  numbers  indicated  by  Eq.  (A-l)  suggests 
that  the  disturbances  will  be  weak  and  linear  theo  ry  may  be  employed.  Let 
us  first  trace  the  acoustic  waves  propagating  from  wall  to  determine  the 
nature  of  the  shock  interaction  zone.  These  waves  are  illustrated  in 
Fig.  A-3.  Wave  A is  the  leading  wave  and  propagates  against  the  mean 
gas  velocity.  Wave  B is  the  trailing  wave  and  propagates  immediately 
behind  the  shock.  The  two  waves  represent  the  limits  of  ihe  zone  of 
influence  of  the  corner. 

To  trace  the  waves  through  the  shock  profile,  we  express  the  wave 
velocity  in  terms  of  the  speed  of  sound  and  the  mean  gas  velocity. 


and 


where 


R 


^2 

cit 


- v Y RT 


(A-  2) 


dR 

l 

dt 


u - V y RT 


(A-3) 


R 


B 


B dt 


Vyrt  , 


(A-4) 


R. 


Rs  (t) 


The  initial  conditions  to  be  imposed  on  Eqs.  (A-2)  through  (A-4)  are:  At 

t - t . 8 8 =9  and  R = R„  = R where  t and  R are  related  by  the 

o’  ABC  ABo  o o 

arrival  of  the  shock  at  R , 

o 


R = At 


2/3 


A-  5 


Equation  (A- 5)  is  illustrated  in  Fig.  A-4.  Wave  A is  initiated  in  a super- 
sonic flow  and  travels  outward,  but  with  a velocity  less  than  the  gas  velocity. 
Eventually  the  subsonic  portion  of  the  shock  profile  overtakes  the  wave  and 
allows  the  wave  to  propagate  to  p = 0.  Since  most  of  the  thrust  is  obtained 
from  the  gas  between  r\-  0.6  and  q = 1 , we  are  interested  in  wave  A between 
T = 1 and  Tj  = 3 only.  The  angular  propagation  of  wave  A is  obtained  from 
Eq.  (A-2).  For  R «=  Rq,  integration  yields 


0 

'A 


V 2y  (V- 1)  / 2/3 

(Y  + l)  \ 1 


(A-6) 


and  § - e is  illustrated  in  Fig.  A-5.  Note  that  in  the  time  of  interest 

C A ^ 

(t  s 3),  wave  A propagates  less  than  30  from  the  wall.  Hence,  for  cone 
angles  greater  than  30°,  wave  A cannot  form  a normal  shock  which  would 
influence  that  portion  of  the  gas  from  which  the  dominant  thrust  is  obtained. 

Before  we  can  verify  the  schematic  of  the  shock  structure  as  illus- 
trated in  Fig.  A-  1 , the  motion  of  the  trailing  wave  must  be  determined. 
Comparing  Eqs . (A-2)  and  (A-4),  we  note  that  deA  a d0B  because 

R aR  . Hence,  shock  A propagates  further  from  the  wall  than  shock  B 
B A 

as  indicated  in  the  schematic.  However,  we  must  still  demonstrate  that 
wave  B can  escape  the  wall.  That  is,  we  must  demonstrate  that  wave  B 
forms  a lamda  shock  rather  than  an  attached  reflection.  Integration  of 
Eq.  (A-4)  yields  G^T^). 


i 


I 


A- 7 


0 S 

B C 


3 (y  + 1) 


(A  - 7) 


If  we  define  0 ^ as  the  angle  v/here  the  initial  shock  intersects  the  wall, 

then  the  simple  geometrical  relationship  illustrated  in  Fig,  A-3  yields 


(R  (t)  - R ) b ..  ? (6  - 6 ) R (t) 

s o wall  C wall  s 


wall  C y T 2/3  J "wall 


(A-  8) 


0B  and  9wall  are  compared  in  Fig.  A-6.  For  6 ^ £ 30  , the  wave  propagates 

faster  than  the  wall  and  must  form  a detached,  or  lamda  shock  as  illustrated 
in  Fig . A- 1 . 


LOSS  IN  THRUST  DUE  TO  SHOCK  REFLECTIONS 


With  the  qualification  that  Fig.  A-l  is  a meaningful  schematic  of  the 

shock  structure  for  0„  > 30°  and  6 < 30°,  we  will  now  proceed  to  analyze 

C wall 

the  effect  of  this  shock  structure  on  the  specific  impulse  of  the  pulsed  laser 
propelled  rocket.  We  shall  calculate  the  velocity  perturbation,  Au»  created 
by  the  shocks  and  demonstrate  that  the  mass  average  fiu  is  very  small  with 
respect  to  the  mass  average  of  u.  Hence,  the  effect  of  the  shock  structure 
is  a small  percentage  loss  in  thrust.  This  will  represent  an  over  estimate 
of  the  losses  because,  once  the  velocity  perturbation  is  created,  it  does 
not  "freeze'1  as  will  be  assumed.  Rather,  the  velocity  perturbation  will 
decrease  as  the  disturbed  region  expands  to  the  nozzle  exit  plane.  This 
effect  is  neglected  in  the  analysis  and  the  results  are  clearly  an  over 
estimate  of  the  PLP  degradation: 

The  mass  average  of  ^ u is  defined  as 

1 

--  / * fe)  ■ 

0 


(A- 10) 


J 


where  M is  the  total  mass  released  per  pulse 


M = 0 u ft  (r  ) (t  - t) 

p s 

and  ft  is  the  solid  angle  of  the  cone 


ft  = 2TT  (1  - cos  S^) 


A-  1 2 


The  system  degradation  due  to  the  reflection  of  the  blast  wave  from 
the  nozzle  wall  will  now  be  assessed.  We  shall  assume  that  the  intersection 
of  shock  B with  the  blast  wave  propagates  with  the  velocity  of  wave  B. 
Therefore,  the  location  of  point  P,  Fig.  A- 1 , is  known  for  all  time  and 
the  mass  affected  by  wave  B is  known  for  all  time.  The  velocity  perturba- 
tion created  by  the  shock  may  be  estimated  in  the  Newtonian  limit 

(£u)  M-)  v e ..  . (A- ll; 

B \ Y 1 / s wa 1 1 

where  5 is  the  maximum  flow  deflection  angle.  The  actual  flow  deflection 

wall 

angle  is  a function  of  Q, 


6 


6 n + 

wall 


3C 


0. 


and  Eq.  (11)  is  an  overestimate  of  the  velocity  perturbation.  Although  this 
may  be  a large  velocity  perturbation,  the  mass  which  is  affected  is  a 
small  fraction  of  the  total  mass  and  the  mass  average  ,\u  will  be  shown 
to  be  small  with  respect  to  u^. 

The  mass  disturbed  by  shock  B is  expressed  as 

2 

dm  - p,  V 2tt  R sin  0 d0  dt, 
p 1 s s 

where  ^ is  the  density  upstream  of  the  shock.  From  (A- 9)  through  A- 11), 
we  obtain 


A-  13 


and  t the  upper  limit  of  the  time  integration,  corresponds  to  either  the 
u 

time  at  which  6 -»  0,  or  the  "break  through"  time,  which  ever  occurs  first. 

It  is  worth  noting  that  we  have  inherently  assumed  that  the  "break  through" 

radius  is  greater  than  R , otherwise  the  shock  reflection  never  forms. 

For  t > t , i.  e.  for  R (tj  > R , Au  is  non  zero, 
bo  s b o 

To  evaluate  Eq.  (A- 13),  we  make  the  approximations  that 
sin  9 * sin  9 , 

and 

de  = ec  - eB  (.)  . 

Integration  of  Eq.  (A- 13)  yields 
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Tlu-  3> 


(A-  14) 


3 + 


1/3 


lu 


where 


4 sin  9 J~2 y (y-1)  (Y-l)  6 .. 

C wal  1 

3 (y  + l)Z  (1  - cos  0 ) I v'R  (t,  )/R 
C c s b o 


and  represents  either  the  break  through  time 


T,  . = (R  (t  )/R  ) 
11  s b o 


3/2 


or  the  time  at  which  h -»  0 


3 (y+1)  6 


12 


= exp 


wa  11 


,2V2y  (Y-  1) 


which  ever  occurs  first.  An  upper  bound  on  Au  may  be  obtained  by  assuming 

that  T and  T occur  simultaneously.  The  maximum  values  of  Au/u  are 
1112  e 

illustrated  in  Fig.  A-7.  For  wall  deflection  angles  above  20  , PLP  degra- 
dation may  be  significant.  However,  this  may  be  reduced  by  choosing  the 
pulse  repetition  frequency  such  that  T ^ and  T are  not  identical.  Thus, 
the  shock  may  " break  through"  before  the  degradation  is  significant.  This 

is  illustrated  in  Fig.  A-8  for  6 = 30°.  If  the  break  through  radius  is 

wall 

less  than  two  R , less  than  10%  loss  will  result.  Should  we  adjust  the 
o 

frequency  such  that  R (t,  ) ■&  R , shock  B never  exists  and  we  need  consider 
n s b o 

only  the  effects  of  shock  A. 
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Rs<»b>/Ro 


Fig.  A-8  Shock  B Degradation  For  a 30° 

wall 
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The  corner  shock.  A,  will  deflect  each  streamline  by  the  amount  6, 


6=6  „ + e - 9 

wall  C 

where  8 is  evaluated  on  each  streamline.  The  velocity  perturbation  asso- 
ciated with  each  streamline  is,  from  linear  theory, 


Au  _ 6 

u = Jm}  - 1 ' (A' 

Since  Eq.  (A- 15)  becomes  unbounded  as  -»  1,  an  upper  limit  on  Au  must 

be  incorporated.  The  maximum  value  of  Au  occurs  for  a plane  normal 
shock. 


(6 u) 


max 


A realistic  composite  of  Eqs.  (A- 15)  and  (A- 16)  is  obtained  by  taking  the 
reciprocal  of  the  sum  of  the  reciprocals. 


(Au)  = k (9,  q)  u = 
A 


(A-  17) 


The  mass  average  Au  is  defined  by  Eq.  (A-9).  For  shock  A,  we  note  that. 


2 

dm  = ou  (2tt  R sin  8 d8)  dt  , 


where  we  have  assumed  that  shock  A remains  at  r : R(  » which  is  in  keeping 
with  the  behavior  of  wave  A.  The  total  mass  M is  that  mass  which  is  in 
the  pulse  created  by  the  blast  wave.  Hence,  Eq.  (A-9)  is  expressed  as 


tu 


u 

e 


4/3  (Y-H) 

VR  (t  )/R  (1-cos  8 )1 
s b o C 2 


k { © , T|)  fh  sin  9 d T d6 

~ 2 
T 1 


(A-  18) 


I'he  choices  of  T.  , and  T are  straight  forward.  We  cannot  create  a 
lx  1 u 

velocity  perturbation  at  9 until  the  wave  gets  there.  Hence,  corresponds 

to  the  time  at  which  wave  A reaches  6, 


(Y+l)  (9C  - 9) 
V2Y  (Y-l) 


The  upper  limit,  T,  corresponds  to  the  time  at  which  M.  -»  1.  At  r = R , 
1 u ° 

q and  T are  related  by 


q = 


(A-  19) 


Hence,  from  Eqs.  (A-l)  and  (A-19), 


T 

lu 


2 

Y (Y-l) 


3/4 
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The  lower  limit  on  9 is  a little  more  complicated.  One  choice  of  9 

min 

corresponds  to  the  value  of  9 at  which  T = T . This  is  denoted  by 
r IX  lu  ' 

9 . , and  corresponds  to  the  minimum  value  of  9 to  which  wave  A can 

min  1 

propagate. 


where 


9 . = 0„  - t , 

min  1 C 


2 - V2y  (V - 1 ) 


A second  choice  for  9 . corresponds  to  the  value  of  9 for  which  the 

min 

deflection  angle  6,  tends  to  zero. 

6 3 8 - 6 

min  2 C wall 


Clearly,  all  values  of  9 must  exceed  both  9 . , and  9 . Hence 

min  1 min  2 


9 . = Greater  of  (9  . ,8  . ) . 

min  min  1 min  2 


The  integrand  of  Eq.  (A- 18)  is  a function  of  T],  T and  9.  However,  at 

r ~ R , T|  and  T are  related  by  Eq.  (A-19).  Therefore,  Eq.  (A-18) 
o 1 

becomes. 


0 T! 

C max 


2/(Y+l) 


u vR  (tT)/R  I (1-cos  0 ) J J 

e sbo2  C 0,  T)  . 


// 


k (8,11)  f fn)  h n)y-  sin  6 dT  d8 


i min 


where 


A"  20 


and 


T\ 

max 


= T -2/3 

1 1 


The  doable  integration  necessary  to  determine  L u/u  has  been  obtained 

e 

numerically.  Results  are  indicated  in  Fig.  A-9.  The  mass  average 

velocity  perturbation  is  only  1°  . Hence,  the  loss  in  thrust  due  to  the 

corner  shock  can  be  no  more  than  1%.  These  results  are  appropriate 

when  the  break  through  radius  is  identical  to  R . For  R (t,  ) > R , the 

o s b o 

effect  of  shock  A decreases  but  degradation  due  to  shock  B increases. 

The  analysis  is  not  appropriate  for  R (t,  ) < R . However  the  effect  of 

s b o 

shock  A would  increase  in  this  situation  because  the  gas  Mach  number 

would  increase  in  the  expansion  from  R (t.  ) to  R . Although  this  effect 

s b o 

has  not  been  assessed,  we  can  approximate  it.  For  isentropic  expansion 

2 

in  a conical  nozzle,  density  decreases  as  1/r  . Hence,  the  speed  of 

sound  decreases  as  r^  . Figure  A-9  indicates  that  the  mass  average 

velocity  can  increase  by  a factor  of  two.  Hence,  for  an  expansion  from 

R (t  ) to  R the  Mach  number  and  velocity  perturbation  can  increase  by 
s b o ^ 

a factor  of  2 (R  /R  (t,  ))  . The  total  degradation  due  to  a 30  w-all  is 

o s b 

illustrated  in  Fig.  A-10.  For  R (t,  )/R  between  0.  1 and  1.0,  less  than 

s b o 

10T  degradation  will  result. 

In  conclusion,  we  have  assessed  the  problem  of  P LP  degradation 
by  blast  wave  reflection  from  the  nozzle  walls.  We  have  shown  that  under 
some  circumstances,  wrall  deflection  angles  of  the  order  of  30°  may  cause 
serious  degradation.  However,  if  we  properly  tune  the  pulse  repetition 
frequency  such  that  the  break  through  radius  is  between  R ^/10  and  R , it 
appears  that  the  degradation  will  be  well  below  10%.  If  R and  R (t ^ ) 
cannot  be  maintained  in  this  ratio,  degradation  may  be  serious,  but  a 
more  refined  model  would  be  necessary  in  order  to  evaluate  the  effect. 
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